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AN  ANALYSIS  OF  GAS  PHASE  ETHANOL/WATER  CHEMISTRY 

FOR  DIAMOND  CVD 

I.  INTRODUCTION 

Diamond  possesses  unique  physical  properties  which  make  it  ideal  for  a  number  of 
technological  applications  including  coatings  for  cutting  tools,  thermal  spreaders  for 
electronic  devices,  high-transmittance  optical  windows  and  high-performance  semi¬ 
conductors.  Previously,  most  artificial  diamond  was  produced  by  high-pressure,  high- 
temperature  techniques,  but  with  the  advent  of  chemical  vapor  deposition  (CVD), 
this  technique  has  become  the  dominant  way  to  grow  high-quality  diamond  films. 
(For  a  more  detailed  treatment  of  diamond  CVD  technology,  see  the  review  article 
by  Celli  and  Butler  [1]).  Early  diamond  CVD  reactors  used  gas  mixtures  such  as 
CO/CO2  [2]  and  CO/H2  [3],  while  later  generations  of  reactors  relied  on  mixtures  of 
hydrogen  and  a  carbon  source,  typically  methane.  In  recent  years  there  has  been  a 
return  to  feedstocks  containing  oxygen. 

An  analysis  of  a  large  number  of  diamond  CVD  experiments  by  Bachmann  tt  al. 
[4]  shows  that  in  nearly  all  reactors  which  were  successful  in  growing  diamond,  the 
compositions  of  the  feed  gases  fell  within  a  well  defined  region  of  a  C/H/0  ternary 
diagram.  Bachmann  et  al.  define  the  following  variables  as  the  coordinates  of  their 
diagram; 


^0=1:  =  [0]/([0]  -h  [^])  (1) 

=  [H]I{[H]  +  [C])  (2) 

^c=i;  =  [C]/(IC]  -k  [0]),  (3) 

and  show  that  the  diamond  growth  region  is  roughly  bracketed  by  the  lines  Xc-x  = 
0.57  and  Xcx  =  0.50.  Their  survey  suggests  that  it  is  not  the  specific  inlet  gases 
themselves  that  are  important,  but  rather  their  overall  C/H/0  ratio. 

The  majority  of  the  CVD  reactors  analyzed  in  ref.  [4]  use  H2  as  an  inlet  gas. 
Although  diamond  films  have  been  grown  quite  successfully  using  H2  as  a  feedstock. 
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one  of  the  serious  drawbacks  is  the  hazard  of  explosion  associated  with  the  storage 
and  handling  of  H2,  as  well  as  of  the  typical  carbon  sources  CH4  and  C2H2.  It  would 
be  advantageous  to  the  CVD  growth  of  diamond  films  to  develop  reactors  which  yield 
high  growth  rates  and  utilize  more  benign  gases.  Recent  experiments  have  demon¬ 
strated  growth  from  methanol  [5],  ethanol/water/hydrogen  [6],  and  EtOH/water  [7]. 

In  this  paper  we  focus  on  the  thermodynamics  and  chemical  kinetics  of 
Et0H/H20  under  conditions  of  constant  pressure  and  temperature.  Previous  stud¬ 
ies  of  C/H/0  systems  include  equilibrium  calculations  of  Et0H/H20  [8]  and  kinetic 
studies  of  CH4/O2/H2  [9-12],  C2H2/O2/H2  [12],  and  CH3OH/H2  [13].  The  calcula¬ 
tions  that  we  perform  are  not  for  any  specific  CVD  reactor  design  or  heating  mecha¬ 
nism,  but  should  be  valid  in  all  cases  where  the  neutral  chemistry  is  dominant.  This 
assumption  is  clearly  not  valid  for  growth  techniques  such  as  electron-cyclotron  reso¬ 
nance  (ECR)  where  the  operating  pressure  is  in  the  mTorr  range  and  e"-impact  is  the 
dominant  H2  dissociation  mechanism,  but  should  be  at  least  partially  valid  for  most 
reactors  operating  at  pressures  greater  than  one  Torr.  Although  it  should  be  possible 
to  grow  diamond  using  a  variety  of  different  alcohols,  we  have  restricted  our  attention 
to  EtOH  since  Et0H/H20  is  the  simplest  alcohol-water  system  where  we  are  able  to 
tune  the  parameter  Xc-.t,  so  as  to  sample  a  region  of  the  ternary  diagram  [4]  from 
“no  growth”,  through  “diamond  growth”,  to  “non-diamond  carbon  growth”.  Using  a 
methanol/water  mixture  it  is  only  possible  to  sample  the  region  of  the  phase  diagram 
below  the  “CO  line”  corresponding  to  Xcx  <  0.5.  An  additional  advantage  of  using 
EtOH  in  our  study  is  the  availability  of  a  reliable  reaction  set  for  EtOH  combustion. 
We  expect  that  the  analysis  of  Et0H/H20  presented  in  this  paper  should  be  easily 
extended  to  other  alcohol-water  systems. 

The  remainder  of  this  paper  is  organized  as  follows.  In  the  next  section  we  de¬ 
scribe  both  the  chemistry  reaction  set  and  the  approach  used  to  perform  the  chemical 
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kinetics  calculations.  The  results  are  presented  in  section  III  with  an  emphasis  on  the 
sharp  divergence  in  species  concentrations  with  inlet  composition,  the  dependence  of 
the  H/H2  ratio  on  operating  conditions,  the  CO/CO2  equilibrium  and  the  ratio  of  Ci 
to  C2  species.  Finally,  in  section  IV  we  summarize  the  results. 

II.  MODEL 

A.  Chemistry  reaction  set 

The  results  of  the  chemical  kinetics  calculations  presented  in  this  paper  were  ob¬ 
tained  using  a  hybrid  chemistry  scheme  constructed  from  the  Miller-Melius  acetylene 
combustion  set  (MM-C2H2)  [14]  and  the  Lawrence  Livermore  National  Laboratory 
EtOH  combustion  set  (LLNL-EtOH)  [15].  Our  hybrid  chemistry  scheme  contains  a 
total  of  249  reversible  elementary  reactions  involving  54  species.  Both  sets  express 
the  forward  reaction  rates  in  the  Arrhenius  form 

k  =  (4) 

and  include  enhanced  third-body  and  pressure  fall-off  reactions.  The  LLNL-EtOH 
set  explicitly  specifies  the  reverse  rate  Arrhenius  coefficients  while  the  MM-C2H2 
set  uses  reverse  reaction  rates  calculated  from  the  equilibrium  constants.  Although 
it  is  generally  considered  unwise  to  tamper  with  reaction  schemes  that  have  been 
optimized  for  the  study  of  a  particular  system,  we  feel  we  have  taken  great  care  in 
combining  and  drawing  on  the  strengths  of  the  two  sets.  The  LLNL-EtOH  scheme 
would  appear  to  be  the  natural  choice  for  this  study  since  it  does  a  very  detailed 
treatment  of  EtOH  chemistry,  but  EtOH  combustion  takes  place  at  temperatures  (T 
<  2000  K)  lower  than  those  used  in  many  CVD  reactors  (T  >  3000  K).  Consequently, 
the  chemistry  of  species  such  as  C  and  C2,  which  are  typically  only  found  at  high 
temperatures,  is  not  considered.  (Thermodynamic  equilibrium  calculations  of  EtOH- 
rich,  Et0H/H20  mixtures  show  that  the  mole  fractions  of  C  and  C2  do  not  exceed 
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1  X  10"^°  for  temperatures  less  than  ~  2000  K  at  0.001  atm.  and  ~  2500  K  at 
1.0  atm.)  The  MM-C2H2  set  on  the  other  hand,  which  adequately  describes  the 
high-temperature  chemistry  of  these  carbon  species,  contains  no  EtOH  chemistry.  In 
combining  the  two  schemes  we  have  taken  all  of  the  elementary  reactions  of  the  MM- 
C2H2  set  and  appended  28  reactions  from  the  LLNL-EtOH  set  to  provide  a  pathway 
from  EtOH  to  species  that  are  treated  in  the  MM-C2H2  set. 

Since  our  CVD  reactor  temperatures  are  considerably  higher  than  those  encoun¬ 
tered  in  typical  combustion  applications,  there  can  be  significant  differences  between 
the  reaction  paths  that  are  important  in  the  two  ca^es.  For  example,  the  unimolec- 
ular  dissociation  reaction  described  below  for  CH3CHO  consumption  is  completely 
dominant  only  at  very  high  temperatures.  As  another  example,  the  H2  dissociation 
reaction  which  is  normally  unimportant  in  combustion  plays  a  crucial  role  under  CVD 
reactor  conditions.  Throughout  the  paper  it  is  important  to  keep  in  mind  the  differ¬ 
ences  between  the  CVD  and  combustion  conditions  and  their  effects  on  determining 
the  dominant  reaction  paths. 

The  LLNL  reaction  scheme  contains  23  reactions  that  directly  involve  EtOH. 
Although  a  rate-of-production  analysis  at  a  range  of  pressures  from  0.001  to  1.0  atm 
and  temperatures  from  2500  K  to  4000  K  indicate  that  only  13  of  the  reactions  make 
any  detectable  contribution  to  the  EtOH  chemistry,  all  23  reactions  were  included 
in  our  reaction  set.  The  only  products  of  the  EtOH  chemistry  that  are  not  present 
in  the  Miller-Melius  set  are  the  species  PC2H4OH,  SC2H4OH,  CH3CHO,  and  C2H5O. 
The  prefixes  p  and  s  on  the  C2H4OH  species  indicate  that  an  H  atom  has  been 
abstracted  from  the  primary  or  secondary  C,  respectively.  In  order  to  merge  the  two 
reaction  schemes,  mechanisms  have  to  be  included  which  break  down  these  products 
into  species  that  are  treated  by  the  MM-C2H2  reaction  set.  Besides  the  reverse 
reactions  which  convert  PC2H4OH  back  into  EtOH,  the  only  other  reaction  involving 
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this  species  is  the  unimolecular  dissociation 


PC2H4OH  =  C2H4  +  OH.  (5) 

Both  reactions  involving  SC2H4OH  lead  to  the  formation  of  CH3CHO, 

SC2H4OH  +  M  =  CH3CHO  +  H  +  M,  (6) 

SC2H4OH  +  02  =  CH3CHO  +  HO2.  (7) 

The  LLNL  set  contains  a  number  of  reactions  involving  CH3CHO  and  C2H5O,  but 
at  the  range  of  operating  conditions  considered,  a  rate-of-production  analysis  shows 
that  the  only  reactions  that  make  any  detectable  contribution  to  their  consumption 
are  the  unimolecular  dissociation  reactions 

CH3CHO  =  CH3  +  HCO,  (8) 

and 

C2H5O  =  CH3  +  CH2O.  (9) 

Both  of  the  reaction  mechanisms  used  to  construct  our  model  were  developed 
originally  for  combustion  applications  at  temperatures  significantly  lower  (T  <  2000 
K)  than  those  encountered  in  typical  CVD  reactors  (T  >  3000  K).  For  the  purpose  of 
the  present  study,  we  have  assumed  that  the  temperature  dependence  of  the  reaction 
rates  in  these  mechanisms  can  be  extended  to  these  higher  temperatures.  In  most 
cases,  higher  temperatures  simplify  reaction  mechanisms  by  reducing  the  number  of 
reaction  paths  which  are  important,  but  details  of  the  reaction  mechanism,  partic¬ 
ularly  those  elements  found  to  be  most  important,  should  be  investigated  at  these 
higher  temperatures  in  order  to  refine  the  model  for  future  detailed  studies. 
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B.  ChemicaJ  kinetics  approach 

All  of  the  chemical  kinetics  calculations  were  performed  using  the  Sandia  Chemkin 
program  Perfectly  Stirred  Reactor  (PSR)  [16,  17],  The  basic  assumption  of  PSR  is 
that  the  mixing  of  reactants  is  so  complete  that  the  conversion  from  reactants  to 
products  is  totally  dominated  by  the  chemical  reaction  rates  rather  than  diffusion, 
convection  or  other  transport  processes.  The  basic  equations  used  in  PSR  are  species 
conservation 


m{Yk-Y:)-UkWkV  =  (10) 

and  energy  conservation 

+  g  =  (11) 

it=i 

where  m  is  the  mass  flow  rate  through  the  reactor,  V  is  the  reactor  volume,  and 
Q  IS  the  reactor  heat  loss.  The  subscripted  variables  Yk,  Wk,  and  t^k  are  the  mass 
fraction,  molecular  weight  and  production  rate,  respectively,  of  species  k.  The  starred 
variables  represent  inlet  conditions  and  the  summation  in  the  energy  equation  is  over 
the  number  of  species  K.  The  mass  flow  rate  through  the  reactor  is  related  to  the 
residence  time  r  in  the  reactor  by  the  equation 


where  p  is  the  mass  density  in  the  reactor.  The  steady-state  conservation  equations 
are  solved  using  a  hybrid  Newton/time-integration  method.  The  PSR  program  can  be 
used  m  one  of  two  modes.  Either  the  temperature  and  pressure  of  the  reactor  can  be 
specified,  in  which  case  the  K  species  conservation  equations  are  solved,  or  the  reactor 
heat  loss  and  pressure  are  specified  and  the  AT  -fl  species  and  energy  equations  are 
solved.  In  all  of  our  simulations  we  assume  that  the  reactor  temperature  and  pressure 
are  known  and  use  PSR  to  calculate  the  species  mole  fractions  in  the  reactor  as  a 
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function  of  time  over  the  range  from  1  x  10"^  s  to  1  x  10^  s.  Since  by  default  PSR 
uses  the  results  of  a  thermodynamic  equilibrium  calculation  as  the  first  guess  for  the 
outlet  gas  composition,  it  is  easy  to  verify  that  the  chemical  kinetics  calculations 
agree  with  thermodynamics  in  the  limit  of  long  residence  times.  For  nearly  all  of  our 
calculations  the  thermodynamic  equilibrium  conditions  are  closely  approached  for  a 
residence  time  of  1000  s. 


III.  RESULTS 

A.  Divergence  of  species  concentrations 

One  of  the  most  striking  features  of  the  high- temperature  EtOH/HjO  chemistry 
is  the  sharp  divergence  in  certain  species  concentrations  as  a  function  of  inlet  gas 
composition.  For  [Et0H]/[H20]  ratios  below  1.0,  corresponding  to  [C]/[0]  <  1,  the 
species  0,  OH,  O2,  H2O,  and  CO2  are  present  in  significant  quantities  (mole  fractions 
>  1  X  10  ®).  As  the  inlet  gas  ratio  is  changed  from  below  to  above  1.0,  the  mole 
fractions  of  these  oxygen-containing  species  drop  dramatically  while  the  mole  frac¬ 
tions  of  the  carbon-bearing  species  C,  CH,  CH2,  C2,  CjH,  and  C2H2,  and  at  higher 
pressures,  CH3  and  CH4,  rise  sharply.  The  mole  fractions  of  the  species  CO,  H,  and 
H2,  in  contrast,  are  relatively  insensitive  to  the  inlet  gas  composition  and  comprise 
the  bulk  of  the  outlet  gases  for  the  full  range  of  conditions  studied.  The  mole  fractions 
of  0,  C,  CO,  H  and  H2  are  shown  in  Fig.  1  at  T=3000  K  and  P=.001  atm.  for  inlet 
composition  ratios  0.5  <  [Et0H]/[H20]  <  1.5. 

The  key  to  understanding  the  divergence  in  species  concentrations  lies  in  the 
extreme  stability  of  CO  [19],  which  guarantees  that  eis  much  CO  as  possible  is  formed 
from  the  available  C  and  0.  For  an  EtOH/H20  ratio  exactly  equal  to  one,  the  C/0 
ratio  is  also  equal  to  one  and  virtually  no  C  or  0  atoms  remain  free  to  form  the 
other  oxygen-  or  carbon-bearing  species  mentioned  above.  At  inlet  composition  ratios 
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greater  than  one,  nearly  all  of  the  0  is  tied  up  in  the  formation  of  CO,  leaving  an 
excess  of  C  atoms  to  form  the  hydrocarbon  and  bare  carbon  species.  Conversely,  at 
inlet  ratios  below  one,  nearly  all  of  the  C  is  tied  up  forming  CO,  leaving  an  excess  of 
0  atoms  to  form  the  oxygen-bearing  species. 

The  strong  dependence  of  species  concentrations  on  inlet  gas  composition  ob¬ 
served  in  Figs.  1  are  a  generic  property  of  the  Et0H/H20  system  and  can  be  pre¬ 
dicted  solely  on  the  basis  of  thermodynamics  [18],  but  chemical  kinetics  calculations 
yield  valuable  information  regarding  the  dependence  of  the  details  of  the  divergence 
on  the  operating  conditions.  To  gain  a  deeper  understanding  of  the  divergence  it  is 
useful  to  look  at  the  time  dependence  of  the  relative  concentration  which 

we  define  as 


(kiri,r2,t)  =  Xkiri,t)/Xk{r2,t).  (13) 

In  this  equation  Xk{ri,t)  and  Xk{r2,t)  are  the  mole  fractions  of  species  k  at  time  t 
and  inlet  Et0H/H20  ratios  ri  and  r2,  respectively.  If  values  for  rj  and  r2  are  chosen 
to  be  slightly  above  and  slightly  below  one,  respectively,  then  significant  deviations 
of  ^2)  t)  from  one  are  indicative  of  the  divergence.  An  advantage  of  studying  (^k 
instead  of  Xk  in  this  context  is  that  it  is  easier  to  compare  time  scales  and  strengths 
of  the  divergences  for  different  species  regardless  of  their  abundances. 

Figures  2  show  ^ic(1.05, 0.95,  f)  for  a  temperature  of  3000  K  and  a  pressure  range 
of  0.001-0.1  atm.  for  all  of  the  major  species  except  H,  H2  and  CO.  Several  points 
regarding  Figs.  2  are  immediately  obvious.  First,  the  divergence  in  concentration  for 
the  different  species  all  occur  at  roughly  the  same  time.  Second,  groups  of  species  can 
be  identified  which  show  nearly  identical  behavior.  These  groups  are  the  C2  group 
(C2,  C2H,  C2H2),  the  Oi  group  (0,  OH,  H2O,  CO2),  and  02-  At  high  pressures  the 
Cl  group  (C,  CH,  CH2,  CH3,  CH4)  is  clearly  evident,  but  with  decreasing  pressure  the 
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c*  for  these  species  show  significant  variation.  Third,  the  C2  and  O2  species  typically 
show  a  stronger  divergence  than  the  Ci  and  Oi  species.  We  classify  CO2  as  an  Oi 
species  for  the  purpose  of  this  study  because  CO,  once  formed,  behaves  essentially  cis 
an  indivisible  unit. 

Because  the  divergence  is  a  gradual  process  and  shows  some  variations  for  different 
species,  it  is  impossible  to  pinpoint  an  exact  time  at  which  the  divergence  occurs. 
Nonetheless,  we  find  it  useful  to  define  an  induction  time,  Tj  ,  which  make  it  possible 
to  identify  trends  in  the  divergence  of  species  concentrations  as  a  function  of  operating 
conditions,  tj  is  defined  somewhat  arbitrarily  as  the  time  at  which  the  Ck  either  exceed 
10.0  or  drop  below  10“T  We  find  that  tj  scales  inversely  with  pressure,  showing  only 
small  deviations  at  higher  pressures,  suggesting  that  the  overall  rate-controlling  steps 
in  our  chemistry  scheme  are  two-body  rather  than  three-body  reactions,  tj  is  shown 
in  Fig.  3  as  a  function  of  pressure  for  several  representative  species. 

In  an  effort  to  identify  the  overall  rate-determining  steps  for  the  divergence,  we 
have  performed  a  sensitivity  analysis  of  the  dependence  of  induction  time  on  the 
Arrhenius  pre-exponential  factors  of  the  elementary  reactions.  At  very  short  residence 
times,  regardless  of  the  inlet  gas  composition  or  reactor  conditions,  large  amounts  of 
CO,  CH3,  CH4,  C2H2,  H2,  and  H2O  are  present.  The  controlling  factor  in  determining 
the  induction  time  should  be  the  rate  of  conversion  of  the  inlet  gases  to  the  principal 
products  CO,  H,  and  H2.  Plots  of  the  species  mole  fractions  as  a  function  of  time 
show  a  correlation  between  the  formation  of  H  atoms  and  the  divergence.  This  finding, 
together  with  the  results  of  the  sensitivity  analysis  shown  in  Table  1,  suggests  that 
the  induction  time  is  most  sensitive  to  the  reaction 

H2-HV1  =  H-|-H-fM.  (14) 

The  induction  time  is  also  found  to  depend,  to  a  much  lesser  degree,  on  the  rates  of 
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the  reactions  +  M  =  C2H  +  H  +  M  and  CO  +  =  HCO  +  H.  Further  analysis 

indicates  that  H  atoms  drive  the  divergence  primarily  through  abstraction  reactions. 
We  have  carried  out  calculations  in  which  we  first  removed  all  reactions  involving  H 
atoms,  with  the  exception  of  the  dissociation/recombination  reactions  listed  above 
and  the  H-atom  abstraction  reactions  involving  EtOH  and  its  derivatives  (PC2H4OH, 
SC2H4OH,  CH3CHO,  and  C2H5O),  and  then  systematically  replaced  reactions  until 
the  divergence  occurred.  Restoring  long-chain  hydrocarbon  reactions,  C„,n  >  3,  and 
reactions  involving  the  breaking  of  C-C  bonds,  be.,  C2H2  -(-  H  =  C  -f  CH3,  to  the 
chemistry  set  has  very  little  effect  on  the  divergence.  It  is  not  until  the  hydrogen 
abstraction  reactions  involving  Ci  and  C2  species  are  replaced  that  the  divergence 
occurs  on  the  timescales  found  using  the  full  chemistry. 

These  results  lead  us  to  believe  that  the  behavior  of  other  alcohol/H20  systems 
should  be  very  similar  to  that  of  EtOH/HiO.  The  rates  for  the  initial  steps  in 
the  dissociation  of  all  alcohols,  the  removal  of  an  H  atom  from  either  a  primary  or 
secondary  carbon,  should  be  very  similar  to  the  H  atom  abstraction  rates  for  EtOH. 
The  additional  reactions  required  for  the  dissociation  of  higher  alcohols,  primarily  the 
breaking  of  C-C  bonds,  would  also  be  expected  to  be  fast  and  not  play  a  critical  role 
in  determination  of  induction  times. 

B.  CO/CO2  equilibrium 

CO2  deserves  some  special  attention  because  it  is  the  only  species,  with  the  ex¬ 
ception  of  short-lived  intermediate  species  and  the  predominating  CO,  which  contains 
both  carbon  and  oxygen  and  is  present  in  significant  quantities  at  high  temperatures. 
The  chemistry  of  CO2  also  is  illustrative  because  it  is  typical  of  those  species  which 
are  formed  in  significant  quantities  initially  and  then  show  a  sharp  increase  or  de¬ 
crease  in  concentration  at  later  times  depending  on  the  inlet  gas  composition.  In 
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the  remainder  of  this  section  we  describe  the  relevant  reactions  which  establish  the 
CO/CO2  equilibrium. 


Various  oxygen-bearing  species  can  attack  CO  to  form  CO2,  but  a  rate-of- 
production  analysis  shows  that  CO2  is  formed  almost  exclusively  by  the  reaction 

C0  +  0H  =  C02  +  H.  (15) 

The  CO2  either  reacts  to  directly  form  CO  or  forms  the  species  HCO  or  CH2O  which 


in  turn  react  to  form  CO  via  the  reactions, 

CO2  +  CH  =  HCO -f  CO,  (16) 

CO2  "h  CH2  =  CH2O  -j-  CO,  (f^) 

C02  +  H  =  HC0  +  0,  (18) 

CO2  "1-^11  =  CO  -h  O  -|-  M,  (19) 

CO2  +  CH2'  =  CH2O  +  CO.  (20) 


In  the  last  reaction  listed  above,  CH2*  is  the  excited  singlet  state  of  CH2.  The  species 
CH2O  is  converted  to  HCO,  primarily  by  interacting  with  H  but  also  by  interacting 


with  OH  and  third  bodies,  through  the  reactions 

CH2O  +  H  =  HCO  +  H2,  (21) 

CH2O  +  OH  =  HCO  -f  H2O,  (22) 

CH2O  -f  M  =  HCO  +  H  +  M.  (23) 

HCO  then  forms  CO  through 

HCO  +  M  =  CO  +  H  +  M,  (24) 

HCO  +  OH  =  CO  +  H2O,  (25) 

HCO  +  H  =  CO  -f  H2.  (26) 
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Through  the  above  sets  of  reactions,  an  equilibrium  will  be  established  between  CO 
and  CO2.  It  is  important  to  note  that  the  only  reaction  pathways  which  exist  for  the 
conversion  of  CO  into  CO2  involve  the  oxygen  containing  species  0,  O2,  OH  and  HO2, 
while  the  conversion  of  CO2  back  to  CO  can  occur  via  reactions  involving  CH,  CH2, 
H,  and  third  bodies.  For  this  reason,  after  long  residence  times  negligible  amounts 
of  CO2  are  found  for  inlet  composition  ratios  above  one,  since  all  of  the  CO2  that  is 
formed  is  rapidly  converted  to  CO. 

C.  H/H2  equilibrium 

Hydrogen  plays  a  critical  role  in  the  growth  of  diamond  films  [19].  Previous 
studies  have  shown  that  both  the  quality  and  the  growth  rates  of  the  diamond  films 
depend  critically  on  the  H  atom  concentration  [19,20].  Although  there  are  a  number 
of  theories  promoting  various  growth  mechanisms,  it  is  almost  universally  accepted 
that  a  crucial  stage  in  the  growth  process  involves  the  abstraction  of  a  hydrogen  atom 
from  the  diamond  surface  through  the  reaction 

CH(S)  +  H  =  C(S,R)  +  H2,  (27) 

where  CH(S)  is  a  hydrogen  terminated  surface  carbon  and  C(S,R)  is  a  surface  radical 
[21].  The  surface  site  created  in  this  way  can  now  accept  one  of  the  carbon  growth 
species.  In  addition  to  initiating  surface  growth,  H  atoms  also  plays  the  key  role  in 
terminating  the  surface  through  the  reaction 

C(S,R)  +  H  =  CH(S).  (28) 

An  analysis  of  the  conditions  which  lead  to  the  growth  of  high-quality  diamond  indi¬ 
cate  that  a  minimum  of  10'‘-24±o-2i  10000-30000)  H  atoms  must  strike  the  diamond 

surface  for  every  carbon  atom  that  is  incorporated  on  the  surface  [22].  Although 
this  paper  is  focusing  primarily  on  gas  phase  chemistry,  it  is  appropriate  to  consider 
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surface  chemistry  from  the  point  of  view  that  these  calculations  could  represent  the 
infinite  limit  boundary  conditions  on  a  ID  transport  or  a  3D  reactive  flow  calculation. 
With  this  emphasis  in  mind,  it  is  now  appropriate  to  consider  what  range  of  residence 
times  and  inlet  gas  compositions  give  rise  to  optimal  H/CHx  boundary  conditions. 
Neglecting  boundary  layer  effects  and  assuming  that  all  CH^  species  have  the  same 
sticking  probability  after  encountering  an  open  surface  site,  we  can  make  the  following 
estimate  of  the  ratio  of  H  atoms  striking  the  surface  to  carbon  atoms  incorporated, 


= 


(29) 


(Ei=o  /  ^rncHs  )'yZ  ’ 
where  7  is  the  sticking  probability  and  Z  is  the  fraction  of  surface  sites  that  are 
open.  We  used  the  values  of  7  =  0.33  and  Z  =  0.05  taken  from  the  work  of  Coltrin 
and  Dandy  [23].  ^  is  displayed  in  Fig.  4  for  a  range  of  inlet  gas  compositions  at 
3.tm.  and  T=3000  K.  At  short  residence  times,  before  the  sharp  divergence 
in  species  concentrations  as  a  function  of  inlet  gas  composition  occurs,  the  ratio  of 
H  atoms  impacting  the  surface  to  carbon  atoms  incorporated  is  much  too  small  for 
high-quality  diamond  growth.  Using  an  Et0H/H20  ratio  of  1.5  yields  values  of  ^ 
that  are  slightly  below  the  critical  value  for  diamond  growth,  while  an  EtOH/HjO 
ratio  equal  to  1.1  results  in  values  of  ^  consistent  with  diamond  growth.  Although 
Fig.  4  shows  that  for  water-rich  inlet  gas  compositions  the  values  of  ^  are  well  above 
the  critical  value,  it  is  unlikely  that  diamond  growth  can  occur  under  these  conditions 
because  of  the  large  concentrations  of  species  such  as  0,  O2,  and  OH  which  have  a 
tendency  to  etch  the  surface  [24]. 

We  have  carried  out  calculations  on  both  pure  hydrogen  and  Et0H/H20  mixtures 
to  determine  the  H/H2  equilibrium.  The  equilibrium  constant  in  the  pure  hydrogen 
system  is  defined  as  /C,  =  [H]/[H2].  For  the  Et0H/H20  mixtures  we  also  define 
an  effective  equilibrium  constant  A"',  which  is  the  [H]/[H2]  ratio  that  is  established 
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when  the  entire  reaction  set  is  considered.  Over  a  large  range  of  pressures,  tempera¬ 
tures,  and  inlet  compositions  the  effective  equilibrium  constant  lies  within  the  range 
l.2lK(.q  <  <  l.SlA'e?'  Thus  the  presence  of  the  other  species  in  the  Et0H/H20 

system  shifts  the  H/H2  equilibrium  slightly  in  the  direction  of  atomic  hydrogen.  Al¬ 
though  the  change  in  the  H/H2  equilibrium  is  small,  it  is  advantageous  that  the  shift 
is  towards  higher  H  concentrations  because  of  the  critical  role  of  H  atoms  in  dia¬ 
mond  growth.  The  relative  insensitivity  of  to  inlet  composition  suggests  that  the 
species  CO,  which  is  itself  insensitive  to  inlet  composition,  is  principally  responsible 
for  the  shift  in  H/H2  equilibrium.  Rate-of-production  analysis  shows  that  CO  drives 
the  equilibrium  in  the  direction  of  excess  H  through  the  reaction  sequence 

CO  -h  H2  =  HCO  -f  H  (30) 

HCO-t-M  =  H  +  CO-hM.  (31) 

The  equilibrium  H/Hj  ratios  for  pure  H/H2  and  Et0H/H20  are  shown  in  Fig. 
5.  These  curves  suggest  a  range  of  possible  operating  conditions  for  a  diamond  CVD 
reactor.  For  example,  even  at  very  low  pressures  it  is  unlikely  that  sufficient  H  atoms 
will  be  produced  at  reactor  temperatures  below  2500  K  for  thermal  plasmas  that  are 
neutral  dominated.  It  should  be  noted  that  the  equilibrium  results  shown  in  Fig.  5 
constitute  upper  limits  on  the  H/H2  ratio.  At  short  residence  times,  the  H/H2  ratio 
for  the  Et0H/H20  mixture  is  very  small.  Since  H  atom  production  follows  primarily 
from  the  dissociation  of  H2,  the  equilibrium  ratios  are  approached  from  below  in  the 
limit  of  long  residence  time.  In  addition,  these  ratios  are  very  temperature-sensitive. 
Since  in  a  diamond  CVD  reactor  the  substrate  temperatures  are  typically  around 
1000  K,  significant  recombination  of  H  atoms  can  occur  as  the  atoms  diffuse  from  the 
hot  interior  of  the  reactor  toward  the  cooler  surface. 
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D.  Relative  concentrations  of  Ci  and  C2  species 

Although  a  gradual  transition  from  diamond  growth  to  graphitic  growth  occurs  as 
the  C/H  ratio  increases  [4],  there  is  uncertainty  over  which  C  species  are  responsible 
for  each  type  of  growth.  Some  studies  suggest  that  the  Ci  species  are  necessary  for 
diamond  growth  while  C2  species  lead  to  the  formation  of  graphite.  Coltrin  and 
Dandy  [23]  support  this  view  by  arguing  on  the  basis  of  thermodynamics  that  the 
surface  adsorption  of  C2H2  is  energetically  unfavorable  as  a  precursor  to  diamond 
but  energetically  favorable  for  graphitic  growth.  Frenklach,  on  the  other  hand,  has 
proposed  a  plausible  surface  reaction  scheme  for  the  growth  of  diamond  from  C2H2 
[25,  26].  Rather  than  investigating  the  validity  of  these  and  other  proposals,  we 
shall  confine  ourselves  to  simply  making  predictions  about  the  relative  abundances 
of  the  Cl  and  C2  species  as  a  function  of  reactor  inlet  gas  composition,  temperature, 
pressure,  and  residence  time. 

The  results  of  simulations  at  temperatures  ranging  from  pressures  ranging  from 
0.001-0.1  atm.  show  that  the  dominant  C  species  at  higher  pressures  and  lower 
temperatures  is  C2H2.  As  the  reactor  temperature  is  increased  and  the  pressure  is 
decreased,  atomic  C  becomes  the  most  abundant  species.  A  crossover  from  C2H2  to 
C  dominance  is  observed  at  a  residence  time  of  approximately  10  s  for  a  temperature 
of  3000  K  and  pressures  of  0.001-0.01  atm.  All  of  these  results  are  illustrated  in 
Figs.  6,  7,  and  8.  In  addition  to  the  temperature  and  pressure  dependencies,  the 
inlet  gas  composition  is  also  seen  to  have  an  affect  on  the  ratio  of  Ci  to  C2  species. 
Figure  9  shows  that  the  ratios  [C]/[C2],  [C]/[C2H],  and  [C]/[C2H2]  all  increase  as  the 
[Et0H]/[H20]  ratio  approaches  one  from  above. 

IV.  SUMMARY 

In  our  study  of  the  high-temperature  chemistry  of  Et0H-H20  we  have  made 
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a  number  of  observations  that  may  have  an  impact  on  the  use  of  these  feedstocks 
for  diamond  CVD  reactors.  The  main  finding  of  this  work  is  the  existence  of  the 
sharp  divergence  in  the  species  mole  fractions  as  the  inlet  Et0H/H20  ratio  is  tuned 
from  below  to  above  one.  Our  analysis  indicates  that  the  divergence  is  due  to  the 
stability  of  CO  and  that  the  primary  rate-limiting  step  determining  the  characteristic 
induction  time  is  the  dissociation  reaction  H2  d-  M  =  H  H  -t-  M.  At  short  residence 
times,  before  the  divergence  occurs,  other  oxygen-  and  carbon-bearing  species  are 
present  in  significant  quantities. 

Although  our  simulations  do  not  explicitly  include  surface  chemistry,  we  can  still 
make  rough  predictions  of  the  operating  conditions  that  will  lead  to  diamond  growth. 
Because  of  the  relatively  low  H  atom  mole  fraction  before  the  divergence,  and  the 
prominent  role  of  atomic  hydrogen  in  promoting  diamond  growth  in  CVD  reactors 
[19],  it  is  likely  that  optimal  diamond  growth  will  only  occur  at  longer  residence 
times.  In  this  limit  the  picture  is  completely  different.  For  Et0H/H20  ratios  less 
than  one,  the  mole  fractions  of  the  potential  growth  species  drop  dramatically  after 
the  divergence  occurs,  making  it  unlikely  that  diamond  growth  will  occur.  Indeed, 
the  availability  of  0,  OH,  and  O2  to  attack  the  diamond  surface  suggests  that  any 
diamond  film  will  be  etched  [24].  For  EtOH/HjO  ratios  greater  than  one,  on  the 
other  hand,  the  growth-species  concentrations  remain  high  after  the  divergence  takes 
place  and  nearly  all  of  the  oxygen  is  tied  up  in  CO.  The  abundance  of  H  atoms  after 
this  split  suggests  a  favorable  environment  for  diamond  film  growth. 

Although  we  have  made  no  predictions  concerning  the  roles  of  different  species  in 
diamond  or  non-diamond  carbon  growth,  we  are  able  to  predict  the  mole  fractions  of 
the  various  species  as  a  function  of  operating  conditions  such  as  the  reactor  residence 
time,  temperature,  pressure  and  inlet  gas  composition.  One  of  the  main  findings 
is  a  crossover  from  C2H2  to  C  as  the  predominant  non-CO  carbon  species  as  the 
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reactor  temperature  is  increased  and  the  pressure  is  decreased.  Another  important 
result  is  that  the  ratios  of  Ci  to  C2  species  increase  as  the  Et0H/H20  inlet  ratio 
approaches  one  from  above.  If  models  based  on  diamond  growth  from  Ci  species 
and  non-diamond  carbon  growth  from  C2  species  are  correct,  then  this  result  would 
be  consistent  with  a  transition  from  diamond  to  non-diamond  growth  as  the  relative 
EtOH  concentration  is  increased. 
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TABLES: 


TABLE  1.  Induction  times  ti  for  the  species  C  and  0  under  operating  conditions 
T=3000  K  and  P=.01  atm.  The  induction  times  listed  in  the  first  table  entry  were 
obtained  using  the  original  unmodified  chemistry.  The  subsequent  induction  times 
were  obtained  for  the  Arrhenius  pre-exponential  factor  of  the  corresponding  reaction 
multiplied  by  ten.  The  effects  of  modifying  all  other  reaction  pre-exponential  factors 
are  small  and  are  not  shown. 


chemistry 

unmodified 

^41 

0.28 

H  -f-  H  -f  M  =  H2  -f-  M 

0.17 

0.13 

H  -|-  H  -(-  H2  =  H2  +  M 

0.16 

0.11 

C2H2 -t- M  =  C2H -f  H -b  M 

0.33 

0.29 

HCO -b  H  =  CO  +  Ho 

0.38 

0.25 

20 


Mole  fraction  O 


t{s) 


(a) 

Fig.  1  —  Time-dependent  mole  fractions  of  the  species  (a)  O,  (b)  C,  (c)  CO,  (d)  H,  and  (e)  H2  for 
T=3000  K  and  P  =  .001  atm.  for  a  range  of  inlet  [EtOHJ/pijO]  ratios  ranging  from  0.5  to  1.5. 
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Mole  fraction  C 


0.001  0.01  0.1  1  10  100  1000 

t(s) 


(b) 

Fig.  1  —  (Continues)  Time-dependent  mole  fractions  of  the  species  (a)  O,  (b)  C,  (c)  CO,  (d)  H,  and  (e) 
H2  for  T=3000  K  and  P  =  .001  atm.  for  a  range  of  inlet  [EtOH]/[H20]  ratios  ranging  from  0.5  to  1.5. 
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0.001  0.01  0.1  1  10  100  1000 

t{s) 


(a) 

Fig.  2  —  Relative  concentrations  vs.  time  for  r,  =  1.05,  r-i  —  0.95  and  a  reactor  temperature  of  3000  K  for  pressure 
equal  to  (a)  0.001  atm.,  (b)  0.01  atm.,  and  (c)  0.1  atm.  is  defined  in  Eq.  13  in  section  Ilia. 
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t(s) 


(b) 

Fig  2  -  (Continues)  Relative  concentrations  f*  vs.  time  for  r,  =  1.05,  rj  =  0.95  and  a  reactor  temperature  of  3000  K  for 
pressure  equal  to  (a)  0.001  atm.,  (b)  0.01  atm.,  and  (c)  0.1  atm.  f*  is  defined  in  Eq.  13  m  section  ffla. 
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(C) 

Fig.  2  —  (Continues)  Relative  concentrations  f,  vs.  time  for  r,  =  1.05,  =  0.95  and  a  reactor  temperature  of  3000  K  for 

pressure  equal  to  (a)  0.001  atm.,  (b)  0.01  atm.,  and  (c)  0.1  atm.  f*  is  defined  in  Eq.  13  in  section  IDa. 
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P  (atm.) 

3  —  Induction  times  t,  as  a  function  of  pressure  at  T=3000  K  for  several  representative  species.  Tj  is  defined  as  the 
time  at  which  Neither  exceeds  10.0  or  drops  below  0.01  for  inlet  composition  ratios  of  r,  =  1.05  and  rj  =  0.95. 


[HMH2] 


Fig.  5  —  [H]/[Hii  equilibrium  ratios  in  pure  hydrogen  (solid  lines)  and  EtOH/Hp  (dashed  lines). 
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Mole  fractions 


Fig.  6  —  (Continues)  Mole  fractions  of  carbon-containing  species  at  T=2500  K  and  inlet  [EtOHJ/PjO] 
ratio  1.1  at  pressure  equal  to  (a)  0.001  atm.,  (b)  0.01  atm.,  and  (c)  0.1  atm. 
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C2H2  C  -A-  C2H  -z-  C2 
CH  -X-  CH2-*-  CH3  -d-  CH4 


Mole  fractions 


t(S) 


(C) 


Fig.  7  —  (Continues)  Mole  fractions  of  carbon-containing  species  at  T=3000  K  and  inlet  [E 
ratio  1.1  at  pressure  equal  to  (a)  0.001  atm.,  (b)  0.01  atm.,  and  (c)  0.1  atm. 
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